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ABSTRACT 

Massive B and Be stars produce X-rays from shocks in high velocity winds with temperatures of a few 
million degrees and maximum X-ray luminosities of ~ 10^* ergs“*. Surprisingly, a sub-group of early 
Be stars exhibits > 20 times hotter X-ray temperatures and > 10 times higher X-ray luminosities than 
normal. This group of Be stars, dubbed y-Cas analogs, contains about 10 known objects. The origin 
of this bizarre behavior has been extensively debated in the past decades. Two mechanisms have 
been put forward, accretion of circumstellar disk matter on to an orbiting white dwa rf, or magnetic 
field interaction between the star and the circumstellar disk (|Smith fc Robins^ll999f) . We show here 
that the X-ray and optical emissions of the prototype of the class, y-Cas, are very well correlated 
on year time scales with no significant time delay. Since the expected migration time from internal 
disk regions that emit most of the optical flux to the orbit of the companion star is of several years, 
the simultaneity of the high energy and optical fluxes variations indicates that X-ray emission arises 
from close to the star. The systematic lack of magnetic field detection reported in recent spectro- 
polarimetric surveys of Be stars is consistent with the absence of strong magnetic wind braking in 
these fast spinning stars but put strong constraints on the possible origin of the magnetic field. We 
propose that in y-Gas the magnetic field emerges from equatorially condensed subsurface convecting 
layers the thickness of which steeply increases with rotation rate and that y-Cas and its analogs are 
the most massive and closest to critical rotation Be stars. 
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1. INTRODUGTION 

y-Gas, the central star of the well known W shaped 
northern sky Cassiopeia constellation is a massive 
and hot BO.SIVe star (M = 16±2M0; Teg «27,000K; 
iStee et all (|2012lB surrounded by a circumstellar disk 
whose matter escaped from the star. The disk radiates 
strong emission lines of hydrogen that together with the 
hot stellar continuum define the spectr al class of Be stars. 
y-Cas was the first identified Be star (jSecchilll866l l. 

Although the optical properties of y-Cas do not dis¬ 
play marked peculiarities, the star stands out from the 
crowd of other Be stars by displaying an X-ray lumi¬ 
nosity of a few 10^^ergs“* (0.3—lOkeV), more than 
one order of magnitude brighter than most B and 
Be stars. The thermal X-ray spectrum is dominated 
by an excee dingly high temperature comp onent with 
fcTf»13keV (iLopes de Oliveira et al.l 1201011 . extending 
up to sa 100 keV (|Shrader et al.l 1201511 . compared to 
kT k: 0.5 keV in normal B and Be stars (iCohen et al.l 
1199711 . In addition, the X-ray light-curve of y-Cas ex¬ 
hibits a unique complex pattern of variability involv¬ 
ing slow modulations on time scales of a few hours 
down to s hot noise lik e even ts as short as a few sec¬ 
onds (s ee iWhite et al.l (1198211 : iLopes de Oliveira et al.l 
(j2010ll : iSmith et al.l ( 201211 arid references therein), y- 
Cas experienced a close to normal B star episode 
in the early 1940’s and since then has more of less 
steadily built a ne w decretion disk and increased its 
optical brightness (iDoazan et al.l Il983t iHenrv fc SmithI 


[2013). Low amplitude quasi-periodic oscillat ions (QPOs) 
with periods in the range of 70±20d (ISmith et al.l 
1200611 dominate the optical variability of the star on 
shorter time scales. Galactic X-ray surveys have re¬ 
vealed a gro wing number (at least 10 so far) of y- 


Cas analogs (ILoues de Oliveira et al.l 

20061: IMotch et al.l 

2007l:lNebot Gomez-Moran et al.l2013 

:IR,auw et al.ll2013l: 

Smith et al.||2015f). 


The origin of this outstanding X-ray emission af¬ 
fecting a small fraction of Be stars remains mysteri¬ 
ous. Two distinct hypotheses have been so far pro¬ 
posed. The discovery that y-Cas is a 203.6 d binary 
with a ft! 1M|^ companion star in a nearly circular or¬ 
bit (iHar manec et al.l 120 001: iMi roshni chenko et al.l 120021 : 
iNemravova et al.ll20i2l : ISmith et al.ll2012ll supported the 
idea arising from some similarities to classical Be/X-ray 
binaries, that accretion of matter onto a compact object, 
most probably a white dwarf, cou l d power the u nusual X- 
ray luminosity ([White et al.lll982l : lHaber]|ll995ll . A hand¬ 
ful of Be stars in the Small Magellanic Cl oud undergo 
nova events at sup er soft X-ray energies ([Sturm et al.1 
HoiaEr et al.ll20T3l . This provides strong support to the 
existence of accreting Be -I- white dwarf binaries, which 
according to binary evolution models should outnumber 
the Be -I- neutron st ar systems by an order of magni¬ 
tude ([Raguzovall2001h . Altern atively, based on corre lated 
optical and X-ray variability (iRobinson et al.l[2002h and 
anti-correlation between UV and X-ray fluxes on time 
scales of several hours, ISmith fc Robinson! ([1999( 1 pro- 
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posed a mechanism in which a magnetic field anchored 
in the rotating Be star interacts with its Keplerian cir- 
cumstellar disk. Reconnection events accelerate parti¬ 
cles that produce X-rays when impacting the disk or 
the photosphere (jRobinson fc SmitM l200r)ll . The origin 
of the putative magnetic field remains mysterious since 
the thick differentially rotating convective layer generat¬ 
ing magnetic structures in solar-type stars is absent in 
such massive objects. 

Here we report the discovery of a strong zero delay cor¬ 
relation between the season averaged X-ray and optical 
light-curves of 7 -Cas obtained over the last 14 years. X- 
ray data were obtained from the All-Sky Monitor (ASM) 
data on-board the Rossi X-ray Time Explorer, (RXTE) 
satellite and from the Monitor of All-sky X-ray Image 
(MAXI) on the International Spaee Station. The Au¬ 
tomated Photometric Telescope (APT) at the Fairborn 
Observatory provided the optical time series. 

2 . X-RAY OBSERVATIONS 

The All-Sky Monitor (|Levine et al.lflOQ^ (ASM) ob¬ 
served 7 -Cas from January 1996 until December 2011. 
We downloaded the “definitive” sum band (1.5 —12keV) 
light-curve accumulated over 90s dwells from the ASM 
RXTE databaseQ. All ASM data acquired in 2011 were 
discarded. During the last year before decommission¬ 
ing a highly erratic behavior and an increase by one 
order of magnitude of the count rate errors prevent 
the use of these late data. On the average, 17 7 -Cas 
dwells are available every day. ASM count rates are cor¬ 
rected for a modeled X-ray spectrum derived from all 
the known sources located in the field of view. As a 
consequence, the count rates recorded by the ASM may 
become negative on occasions. For a comparatively faint 
X-ray source such as 7 -Cas about 37% of the sum band 
d wells rates are neg a tive. Following the procedu re used 
in [Trowbridge et ^ (|2007ll and recommended by iLevind 
(j2013[l we included negative count values in our analysis. 
We rebinned the 90 s averages to longer time bins us¬ 
ing a weighted mean process encoded in an ESO-MIDAS 
procedure. 

The MAXI instrument (iMatsuoka et al.l 1200911 is an 
all-sky imager operating in the 2—20keV energy band 
installed on the Japanese Experiment Module of the In¬ 
ternational Space Station (ISS). The sky coverage is of 
85% per 92-minute orbital period and of 95% per day. 
7 -Cas was observed from August 2009 on. We down¬ 
loaded data averaged over one orbital period from the 
MAXI databas^l- Time series counts and associated 
errors were collected for the total 2—20keV band and 
for the 2—4keV, 4—lOkeV and 10—20keV bands. The 
best signal to noise ratio is achieved in the 2—20keV 
band. We rebinned the one orbit data to longer time 
bins using the same weighted mean process as for ASM 
data. MAXI light-curves contain occasional poor qual¬ 
ity data pointfl in the form of flare-like events due to 
an unexpected background increase or dip-like profiles 
caused by shadowing by solar-battery paddles or ISS 

^ http://xte.mit.edu/asmlc/ 

^ http://maxi.riken.jp/:op/. Our analysis is based on data 
downloaded on '2014 May 7. Consequently, these data do not suffer 
from the incorrect storage of the flux errors reported by the MAXI 
team. 

® see light-curve product notes at http://maxi.riken.jp/ 



Figure 1. The relation between the ASM (1.5—12.0 keV) X-ray 
flux and the optical V band flux expected from the decretion disk 
(see text) in the time interval 2003-2011. These data show seasonal 
(yearly) averages, as discussed in the text. The dotted line shows 
the best fit linear relation. Fluxes are given in ergcm“^ s“^. Error 
bars show statistical uncertainties at the one cr level. 

servicing spacecraft docking. In the 5 and 10 d aver¬ 
ages light-curves, we retained all data points with as¬ 
sociated Icr errors below O.OIcts/s and count rates in 
the range of 0.01 to 0.05. We converted ASM and MAXI 
count rates into X-ray fluxes using pimms version 4.7 and 
assuming a spec tral shape similar to that obs erved by 
XMM-Newton in fLopes de Oliveira et al1(|2010D . namely 
Nh = 2 X 10^^cm“^, abundance 0.1 solar and thermal 
plasma emission with kT = 12 keV. We found that a one 
ASM cts/s corresponds to 3.47 x 10“^° ergcm“^s“^ in 
the 1.5 —12.0 keV band while one MAXI broad band 
cts/s/cm^ corresponds to 7.82 x 10“® ergcm“^s“^ over 
the same energy band. 

The combined high sensitivity and large field of view 
offered by recent all-sky X-ray monitors has yielded a 
wealth of high quality information on the variability of 
cosmic X-ray sources on time scales ranging from sev¬ 
eral years down to as short as a few seconds. However, 
their high “grasp” comes at the expense of data reduc¬ 
tion simplicity and on occasions, can be at the origin 
of spurious features. For instance, coded mask instru¬ 
ments such as the ASM are subject to background mod¬ 
eling uncertainties due to the varying merged contribu¬ 
tion of X-ray sources in the field of view. Likewise, slit 
-I- collimator instruments such as the MAXI SSC and 
GSC cameras may also display background and signal to 
noise variatio ns due to the orbital precession of the ISS 
(|Mukail[2015D . We show in the appendix that the ASM 
and MAXI light curves of 7 -Cas display periodicities of 
~ 200 d and ~ 70 d, respectively, and that these periodic¬ 
ities are intrinsic to the all-sky monitors. These two pe¬ 
riods are very unfortunately close to the orbital period 
of 7 -Cas and to the mean period of the quasi-periodic 
oscillations seen in the optical photometry. We can nev¬ 
ertheless conclude, based on the ASM light curve, that 
70 d X-ray quasi-periodic oscillations are likely present 
and, based on the MAXI observations, that the X-ray 
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emission is not modulated at the 203.6 d orbital period 
of 7 -C as, in agreement with the results of iShrader et al.l 

ll20m . 

3. OPTICAL OBSERVATIONS 

The APT optical monitoring consists of 16 observing 
seasons covering the time interval 1997 September 27 
until 2013 October 2. Data are availab l e on-l ine from 
iHenrv fc SmithI (|2012f) and from iHenrvI (I2014I1 . How¬ 
ever, due to a calibration uncertainty, the first season 
cannot be used for the purpose of long term monitoring. 
The Johnson B and V filter observations of 7 -Cas were 
made with the T3 0.4m APT at Fairborn Observatory in 
southern Arizona. F ull details conce r ning the observa¬ 
tions can be found in iHenrv fc SmithI (| 20 I 2 f) and papers 
listed therein. Briefly, the APT was programmed to ob¬ 
serve 7 -Cas once every two hours during every clear night 
of the observing season (consisting roughly of the month 
of June, interrupted by the Arizona summer monsoon, 
and September-February). In addition, several nights 
were targeted for concentrated observations. In these 
cases the observing sequence in both filters (with a check 
star and comparison star to guard against systematic er¬ 
rors) occurred at a cadence of 8 minutes. The rms errors 
are ± (0.003 — 0.004) mag. The yearly-mean check minus 
comparison differential magnitudes scatter around their 
overall mean with an rms of 0.0019 mag. 

4. X-RAY/OPTICAL CORRELATIONS 

Optical emission lines beca me extremely weak or 
even absent in the early 40’s (jCowlev & MarlboroughI 
Il968f) . At around the same time, 7 -Cas reached the 
faintest V magnitu de historically recorded at V^2.8 
(jPoazan et al.ll98^ . It is thus likely that the star had al¬ 
most completely lost its circumstellar material circa 1942 
and that the V magnitude recorded at that time repre¬ 
sents the brightness of the naked star. Over the time 
interval 1998—2013 covered by APT monitoring, the av¬ 
erage optical emission of the circumstellar disk had in¬ 
creased the brightness of 7 -Cas by about 0.6 magnitude 
in V with respect to naked stellar emission. Accordingly, 
we estimated disk optical emission by subtracting from 
the APT V flux that observed from the naked star. V 
band fluxes were computed assuming an effective filter 
width of 880 A, and a V zero point of 3.75 x 10 “® erg 
s“^ cm“^ (jAllenll20(Hll . The V band disk flux is then 
given by: 

= 3.7510"^ X 880 x (10"*-10"i^ ) ergs“^ cm'^ 

.W 

APT optical data are not evenly distributed in time 
due to the visibility of the star and to the monsoon 
season. In order to build long term optical and X- 
ray light curves, we grouped optical data inside specific 
windows that for each year exclude a few (about 8 %) 
scattered V and B magnitudes obtained before or after 
the main observing season which extends from approxi¬ 
mately September to January. We computed ASM and 
MAXI count rate averages over these time windows so as 
to ensure the best simultaneity between X-ray and op¬ 
tical data. Over the time interval from 1998 until 2012, 
the mean duration of the yearly time bin is 0.43 yr. How¬ 
ever, because of the much coarser coverage obtained dur¬ 
ing the last optical runs, the last data point is averaged 



2000 2005 2010 2015 

Year 

Figure 2. ASM, MAXI and V band from APT monitoring yearly 
season light-curves. Black squares: ASM (1.5 — 12.0 keV) X-ray 
flux. Blue diamonds: MAXI fluxes in the ASM energy range. The 
spurious 200 d ASM modulation has been removed. Red squares: 
V band flux from the decretion disk scaled down to the X-ray flux 
range using the relation derived from ASM and V band data af¬ 
ter 2002 (see text). MAXI data follow the same trend as earlier 
ASM data after 2010. The larger MAXI flux error bars account 
for the 10% cross calibration uncertainty, however, their intrinsic 
mean error is of only 7 x 10“^^ ergcm~^s“^. The origin of the 
discrepancy between ASM and MAXI fluxes in 2009 is unknown 
and since the ASM and MAXI fluxes are identical in 2010, we sus¬ 
pect a calibration problem during the first MAXI operation year. 
Error bars show statistical uncertainties at the one cr level. Errors 
on scaled down V flux are too small (~ 10“^^ ergcm”^ to be 
visible. 



Figure 3. The crosscorrelation function of the ASM and APT 
light-curves shown on Fig. [2] The gray parts correspond to the 
superposition of 50,000 crosscorrelations which were constructed 
from 50,000 pairs of fake light curves. The crosscorrelation of the 
original light curves is shown in black. 


over a 1.34 yr long time interval from May 2012 until Oc¬ 
tober 2013. Fig. [T] shows that a well defined correlation 
between X-ray and disk optical emission exists on time 
scales of 9 years from 2002 until 2010. The relation be¬ 
tween the ASM (1.5 —12.0keV) flux and the optical V 
band emission from the decretion disk computed using 
CD can be simply fitted by a linear relation. Taking into 
account errors on optical and X-ray flux we find: 

^ ^29 (±14) X ±1.75 (±0.02) 10"^ ergs"^ cm” 

( 2 ) 

Extrapolating this relation at null X-ray flux and assum¬ 
ing that a linear regime holds over the whole flux range 
implies a “basal” V band flux equivalent to 85 (± 6 )% of 
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Figure 4. ASM and V band flux light-curves. Times are in years. Black squares: ASM data (1.5 — 12keV) averaged over 10 d time bins 
after removal of the spurious 200 d modulation (see appendix). Red squares represent the 1 d averages of the variable part of the optical 
circumstellar disk emission scaled down to the X-ray flux range (see text). Error bars show statistical uncertainties at the one cr level. Blue 
triangles and associate d error bars mark the full range of X-ray flux detected in 2000 by the PCA on RXTE as based on 30 minute averages 
JRobinson et al1l2002l) . 
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Figure 5. X-ray and optical periodograms. The periodogram 
of the optical V band data (orange contours showing normalized 
powers 15,20 25) is overplotted on the ASM (1.5 — 12.0 keV) peri¬ 

odograms (normalized power range from 0 to 2). All periodograms 
are computed on a one year long interval sliding by 0.1 year and 
Gaussian filtered in the time direction with a sigma of 0.3 year. 
ASM data have a 5 day time resolution while optical measure¬ 
ments are averaged over one day. ASM data have been cleaned 
from the spurious 200 d artifact. Some excess power is detected in 
the period range from 50 d to 100 d. The strongest optical ^80d 
oscillations are weakly correlated with the largest X-ray power. 


the disk V band flux observed in late 2010. Here we as¬ 
sume an error of ± 0.05 mag on the V magnitude of the 
naked star. This suggests that only a relatively small 
fraction of the regions of the decretion disk contribut¬ 
ing to V band emission is related to the production of 
X-rays. Our AFyratio of ss 130 is in agree¬ 
ment with the value of ~ 80 found by iSmith et al.l (j 200 fi[ i 
when correcting for stellar contribution. We then used 
this relation to scale the optical flux of the decretion 
disk to the X-ray range and plot on the same frame the 
ASM, MAXI and APT season averages shown in Fig. ^ 
The light curves confirm the amazingly good correlation 
between X-ray and optical data acquired after 2002, ex¬ 
tending to the time of the MAXI observations. A similar 
correlation probably occurred before 2002 , albeit with a 
slightly shifted zero point. Brighter optical disk emis¬ 
sion is invariably accompanied by a similar increase of 
the X-ray flux. 

We estimated the maximum lag between optical and 
ASM flux variations using a Monte Carlo method sim¬ 
ulating a total of 50,000 different light-curves with the 
same statistical properties as the observational ones. The 
crosscorrelation function of the ASM and APT light 
curves and the results of the Monte Carlo simulations 
are shown in Fig. [3l We find that any delay between 
X-ray and V band flux variations is less than about one 
month at the 3 a level. 

The binned ASM X-ray light-curves suggest the pres¬ 
ence of oscillations on time scales of ^ 50 to 100 d (see 



Figure 6. Power spectra of strips of ASM data. Black: 2000.8 - 
2001.8, green: 2001 - 2001.8, red: 2005.2 - 2006.2 and blue; 2008.5 
- 2009.2. QPOs with periodicities in the range of 70 d to 110 d are 
detected in these particular time intervals. However, QPO periods 
down to about 50 d are also seen on occasion (see Fig. 0. 


Fig. U a,e for instance and Fig. [S]). Similar time scales 
have been seen in the op tical time series collected since 
1996 (|Smith et al.l 1200611 : see also Figs. |4] and [5l Over¬ 
all, the statistical significance of the X-ray QPOs is not 
very high. However, X-ray QPOs exhibit stronger ampli¬ 
tudes during some time intervals. The four power spec¬ 
tra shown in Fig. [S] have bin probabilities < 2% to be 
due to a white noise random fluctuation at the optical 
QPO periods. It seems therefore very likely that the 
X-ray light curve of 7 -Cas is also modulated with quasi- 
periodic oscillations at a mean period very close to that 
reg ularly detected in the o ptical photometry. For exam¬ 
ple, [RoHniQKAE^ l|2002D reported evidence of a ~ 70d 
X-ray modulation in a series of 6 snapshot RXTE/PCA 
observations obtained in 2000. Accordingly, we over¬ 
plot on Fig. m b & c the scaled observed PCA flux on 
the top of the equivalent ASM 10 d broad band aver¬ 
ages. Error bars represent the total range of PCA count 
rates averaged over 30 mn time intervals observed during 
each of these r:; 27h contin uous RXTE/PCA observa¬ 
tions (|Robinson et al.l I2002L Fig. 6). Mid PCA fluxes 
match very well those derived from ASM data provid¬ 
ing an independent confirmation of the reliability of the 
ASM light curves. Although in 2000, the 70 d QPOs are 
detected at both X-ray and optical energies, with hardly 
any time delay, it can be seen that the relatively few 
RXTE/PCA short observations cannot accurately sam¬ 
ple the full complexity of the X-ray light curve. 

In spite of the very likely existence of a Ri 70 d sig¬ 
nal at both X-ray and optical wavelengths, the over¬ 
all crosscorrelation between ASM X-ray fluxes and op¬ 
tical fluxes does not show any significant peak. This 
probably reflects the complex and changing relation pat¬ 
terns between optical and X-ray flux variations on weekly 
timescales. However, on several occasions strong and well 
correlated trains of optical and X-ray modulations are 
clearly detected with small absolute time lags (see e.g. 
Fig. Sic, d, e & i). 

Importantly, the relative amplitudes of the X-ray and 
optical modulations on time scales of a few weeks are the 
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same as those derived from yearly averages, thus hinting 
at a common mechanism for both the long and short time 
scales variations. 

5. DISCUSSION AND CONCLUSIONS 

The optical brightness of 7 -Cas (U ~ 2.4) has turned 
the star into an obvious target for interferometric cam¬ 
paigns. Observations carried out in the visible (Ha con¬ 
tinuum) and ne ar infrared (H b and) between 2010 Aug. 
5 and Nov. 10 (iStee et al.ll2012tl give a disk diameter of 
~ 1.9 D* indicating that most of the visible and near-IR 
disk emission arises from the 0.45 R* wide inner part of 
the disk above the photosphere. This is consistent with 
viscous decretion di sks modeling (|Carciofi fc BiorkmanI 
l2006t iCarciofil 1201111 which indeed predict that for sys¬ 
tems similar to 7 -Cas, V band emission is confined to 
within 1 R* from the B star. Therefore, the « 70 d V 
band QPOs as well as the long term changes must orig¬ 
inate from relatively close to the star. The mechanism 
giving rise to the 70 d QPOs is not known, although the 
possibility that these oscillations are due to a dynamo 
operating in the decr etion disk has been proposed by 
iRobinson et al.l (l2002f l. At any rate, the 70 d periodicity 
is probably not directly related to Keplerian rotation. 
Since 7 -Cas is rotating at or n early at the critical equa- 
torial velocity d^ee et al.ll 2012 tl with a period of « 1.22 d 
(jHenrv fc Smithll2012ll . a 70 d Keplerian rotation period 
is reached at R = (70/1.22)^/^R* rs 15R*, well beyond 
the end of the V band emitting region. The proxim¬ 
ity of the ratio between orbital and mean QPO peri¬ 
ods to 3:1 suggests a possible resonance effect. Indeed, 
low eccentricity systems such as 7 -Cas are expected to 
display moderate disk truncation a t the 3:1 resonance 
radius ([Okazaki fc NegueruelallWiTI l. Models of tidally 
perturbed disks show density spiral patterns that may 
extend to the inner parts of the disk. Such an orbitally 
“locked” mechanism would explain why the 70 d QPOs 
clock is at work ind ependently of disk optical brightness 
(jSmith et al.l 1200^ . Obseryational eyidence in fayor of 
tidally induced effects remain rather scarce. One ex¬ 
ample is the binary Be star tt Aqr, the Ha profile of 
which shows orbital phase-locked Violet/Red lobe yari- 
ations apparently caused by a di sk region of enhance d 
density facing the companion star (IZharikoy et al.l[2013f) . 

Our analysis reveals for the first time the excellent long 
term correlation between the yearly averages of the V 
band brightness of the circumstellar disk and X-ray emis¬ 
sion. X-ray and optical correlated behaviors o f the 70 d 
QPOs were suspected bv IRobinson et all (|2002ll based on 
a coarse X-ray monitoring. We show here that although 
such a strong correlation may happen on occasion, the 
actual correlation pattern is complex. 

The high value of the ratio between optical and X-ray 
flux variations on time scales of years or on time scales 
of the order of the 70 d QPOs, (AFy ss 130) 

rules out a mechanism in which optical light changes are 
due to X-r ay heating. A sim ilar conclusion has been 
reached by iSmith et al.l (119981 1. In principle this does 
not exclude the possibility that a mass ejection eyent 
first detected as a surge in the V band from the inner 
disk migrates to the orbit of the companion star, is ac¬ 
creted and giyes rise to an increase of X-ray luminosity. 
Likewise, the density modulations in the inner part of 



10000 15000 20000 25000 50000 
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Figure 7. The distribution of Be stars in the versus Vrot/Vc 
diagram for stars having inclinations measured by long base inter¬ 
ferometric observations. Green symbol: HD 110432. Red symbol: 
y-Cas. 


the decretion disk seen as 70 d optical QPOs could per¬ 
haps be transported to a putative accreting companion 
and modulate the accretion rate on similar time scales. 
However, the near simultaneity of the long term optical 
and X-ray variations rules out this possibility. 

Viscous disks response t o sustained ni a ss inj ection 
has been inv e stigat ed by I Okazaki et al.l (1200211 and 
iCarciofi et al.l (1201211 amon g others. Dynamical mod- 
els bv iHaubois et al l (1201211 show that the inner parts 
of the disk are first filled in and that at radii equiva¬ 
lent to that of the prim ary Roche lobe (i? s»21R* for 
7 -Cas (IGies et al.ll2007ll l. steady state disk densities are 
only reached after typically 5 to 20 years depending on 
the assumed a viscosity. The dissipation of the decre¬ 
tion disk of the Be star 28 CM a is well described by 
these models (|Carciofi et al.l[20T^ . Unfortunately, to our 
knowledge, only few Be/X-ray systems have long enough 
X-ray and optical light curves to test the validity of the 
dynamical models at relatively large distances of the B 
star. The first example is SXP18.3, a Be/X-ray neu¬ 
tron star binary in the SMC with a 17.79 d orbital pe¬ 
riod which experienced a period of strong X-ray activity 
about 4 years after the be ginning of a major m ass ejec¬ 
tion event seen optically (iSchurch et al.1 [200^ . Nearly 
the same delay between major optical and X-ray out¬ 
bursts may ha ve been seen in the low Lx persistent X-ray 
pulsar X Per (|Li et al.ll2014ll . With an orbital period of 

250 d, an eccentricity of e = 0.11 and an O9.5HIe-B0Ve 
primary, X Per resembles 7 -Cas. This system does not 
exhibit strong X-ray flux variations with periastron pas¬ 
sage. However, no significant correlation with time lags 
lower than a fe w years is seen betwee n optical and X- 
ray light curves (iLutovinov et al.ir2012ll , at variance with 
the observed behaviour of 7 -Cas. Finally, the Be/X-ray 
binary A 0535-1-26 terminated in 2005 its ten year long in¬ 
activity period by a series of giant and normal outbursts 
starting about hve years after the onset of a V band flux 
increase marking the rebuilding of a new dec retion disk 
(|Yan et al.l 120121 : iCamero-Arranz et al.ll2012D . In these 
systems, the long time delays between outbursts in X- 
rays, due to the sudden availability of matter for accre- 
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Table 1 

X-ray (ROSAT) and interferometric rotational parameters of the Be stars plotted in Fig. [3 


HD 

number 

Log(Lx) 
(0.1-2.4 keV) 

Vrot/Vc 

Spectral 

type 

TeS 

(K) 


i 

(deg) 

Ref;^ 

5394^ 

32.13 

1.15 

± 

0.15 

B0.5IVe 

26431 

± 

618 

47 

1 

10516 

< 29.78 

0.82 

± 

0.06 

B1.5Ve 

25556 

± 

659 

72 

1 

50013 

29.85 

0.80 

± 

0.31 

B1.5V 

24627 

± 

590 

31 

2 

120324 

30.45 

0.72 

± 

0.20 

B2IV 

22554 

± 

661 

18 

2 

110432 

31.08^ 

1.00 

± 

0.20 

BO.5-1 

20324 

± 

344 

63 

1 

37202 

< 29.27 

0.71 

± 

0.02 

B2IVp 

19310 

± 

550 

78 

1 

158427 

28.57 

0.90 

± 

0.17 

B3IV 

18044 

± 

310 

39 

2 

142983 

< 30.02 

0.97 

± 

0.12 

B8Ia/Iab 

17645 

± 

554 

57 

1 

91465 

< 29.15 

0.76 

± 

0.15 

B4V 

17389 

± 

415 

45 

1 

25940 

< 29.36 

0.84 

± 

0.22 

B3Ve 

16158 

± 

582 

43 

1 

22192 

< 29.27 

0.78 

± 

0.18 

B5Ve 

15767 

± 

509 

72 

1 

4180 

< 29.83 

0.75 

± 

0.09 

B5IIe 

14438 

± 

272 

57 

1 

37795 

27.75 

0.95 

± 

0.23 

B7IV 

12963 

± 

203 

33 

2 

209409 

< 29.49 

0.74 

± 

0.05 

B7VIe 

12942 

± 

402 

78 

1 

209409 

< 29.49 

0.77 

± 

0.21 

B7IV 

12492 

± 

402 

45 

1 

23630 

< 28.98 

0.79 

± 

0.09 

B7IIIe 

12258 

± 

505 

44 

1 

58715 

< 28.49 

0.86 

± 

0.25 

B8Ve 

11772 

± 

344 

45 

1 

89080 

< 28.87 

0.84 

± 

0.16 

B8III 

11720 

± 

431 

55 

1 


® References, ('ll IBerghoefer et al.l Ill99fili : f2l lCohen et al.l 1119971 ! 

HD 5394 = 7 -Cas 

Lx = 4.2-5.2 X 10®^ergs“l (0.2-12keV) IILopes de Oliveira et al.|[2007l1 


tion at the neutron star orbital distance, and in optical 
signatures of a sustained increase of the decretion density 
close to the star, are consistent with the predictions of 
the dynamical models. Consequently, the maximum one 
month lag between X-ray and optical long term changes 
completely rules out accretion on a compact object in a 
203.6 d orbit as the X-ray emission mechanism in y-Cas. 
Optical variations can be easily understood as changes 
in the emission measure in regions within the inner disk. 
These arguments are consistent with magnetic disk-star 
interactions as the most natu ral X-ray emitting mecha - 
nism, as originally proposed (jSmith fc B.obinsonl 11999(1 . 
We develop this point next. 

The “Magnetism i n Massive Stars” (MiMes) project 
(|Grunhut et al.l[2M^ showed that strong and structured 
magnetic fields of likely fossil origin can be found in 
about 7% of all massive O and B stars. These fossil 
helds may have their origin in a convecting phase oc¬ 
curring during pre-main sequence stage, in an early bi¬ 
nary nierger_ 2 t^ej_ 0 £^be remnant of a primordial field 
(see lWalder et al1l2012L for a recent review). However, 
one of the most intriguing result of MiMes is the to- 
tal lack of magnetic field detection in classical Be stars 
(|Grunhut et al.l[2012ll with upper limits of B ^ 250 G and 
B ^ 500 G for organized and s mall scale fields with filling 
factors of 50% respectively (|Rivinius et al.l [201^ . Im¬ 
portantly, y-Gas also fails to show e vidence of g lobal 
magnetism with similar upper limits (|Wadel [20141) . al¬ 
though the detection of the stellar rotation period in the 
optical photometry (jHenrv fc SmithI 120121) suggests the 
presence of some magnetic spots at the surface of the B 
star. A magnetized wind can efficiently carry away angu¬ 
lar momentum and spin down ste llar rotation with time 
scales of the order of a few Myr (|Petit et al.ll2013ll . sig¬ 
nificantly shor ter than the estim ated ages of y-Gas and 
of its analogs (jSmith et al.ll2012ll . The apparent absence 
of strong and structured magnetic field in Be stars is 
therefore consistent with their high rotation rate. 


Thin subsurface convecting layers may exist in all mas¬ 
sive stars as a result of t he iron peak in stellar opaci¬ 
ties (|Gantiello et al.ll2009tl . The thickness and mass con¬ 
tent of the convective layers increases very substantially 
and are more equatorially co ndensed for stars ro tating 
close to the critical velocity (|Maeder et al.l l2008ll such 
as y-Gas. This mechanism could explain the presence 
of strong, localized magnetic fields connecting the stellar 
equator to the decretion disk and generating the anoma¬ 
lous X-ray emission. Being mostly located in the equato¬ 
rial plane such random magnetic fields may remain con¬ 
cealed from optical spectro-polarimetry observations. 

Long baseline interferometric survey observations of 
Be sta rs h ave been recently repo rted bv iMeilland et all 
(|2012f) and iTouhami et ^ (120131 1. Disk axial ratios pro¬ 
vide information on inclination and together with ob¬ 
served velocities, suggest that on the average, Be stars 
rotate at about 0.8 times the criti cal velocity. S imilar 
conclusions have been reached by iFremat et al.l (j2005f l 
based on the modeling of optical spectra taking into ac¬ 
count the flattening and varying gravitational darkening 
caused by the fast rotation. Interestingly, long base¬ 
line optical interferometric observations of the the two 
brightest y-Gas analogs, y-Gas itself (BO.5 IV-Ve) and 
HD 110432 (BO.5 Hie) confirm that they are both nearly 
critical rotators (Vrnt./ yn = 1.15 ±0.15 and Vrnt,/Vn = 
1.0 ±0.20 respectively (iStee et al.ll2012L [201311 1. We note 
that for thes e two stars, the inclinations computed by 
IFremat et al.l (l2005f ) are substantially higher than those 
derived from long base interferometric observations, sug¬ 
gesting a possible breakdown of the optical spectral mod¬ 
eling at close to critical velocities. 

We show in Fig. [7] the distribution of 18 Be stars with 
interferometrically measured inclinations in the Teff ver¬ 
sus Vr nt./Vn diagram. We e xcluded the ill-defined S Gen 
value ifM eilla,nd et al.l [21)1^ . All Teg were taken from 
IFremat et al.l (j2005| l. All Be stars marked in Fig. [7] have 
been either detected in X-rays or have signihcant upper 
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limits on their X-ray luminosity, either from the ROSAT 
all-sky survey (jBerghoefer et al.il996l : lVoges et al.lll999li . 
or from the ROSA T survey of Be stars conducted by 
iCohen et al.l ()1997[ i. The corresponding values are listed 
in Table [TJ 

Although error bars remain large, Fig. [3 shows that 
both 7 -Cas and its analog HD 110432 are located in a re¬ 
gion of extreme rotation rates and hot effective temper¬ 
atures apparently void of normal Be stars. In particular, 
none of the four highest (> 19,000 K) Be stars with 
Vrot/Vc < 0.9 (HD 10516 = </> Per, HD 50013 = k CMa, 
HD 120324 = /i Cen and HD 37202 = ^ Tan ) exhibit an 
X-ray luminosity in the typical range of 7 -Cas analogfl 

We therefore propose that 7 -Cas-like objects could 
represent the most rapidly rotating of the most massive 
Be stars in which magnetic field generation is the di¬ 
rect consequence of critical rotation. The fa ct that some 
of the 7 -Cas analogs a re blue stragglers (iMarco et al.l 
120091: iSmith et al.ll2012[ l which could have been spun up 
to highest rates through binary evolution provides fur¬ 
ther support to this scenario. However, the very narrow 
range of B0e-B1.5e spe ctral types in which the 7 -Cas phe¬ 
nomenon is observed (|Smith et al.ll201^ certainly calls 
for an explanation. The hot boundary is likely related to 
the disappearance of circumstellar disks at this spectral 
type due to radiative pressure or due to a different evolu¬ 
tion o f the angular momentum (see e.g. lNegueruela et al.l 
I2004L and references therein). The lower mass limit may 
reflect the steep decrease of the size of the convective lay¬ 
ers around B2 and may be related to the observed sharp 
drop of the Lx /Lhni ratio at about the same spectral type 
dCoLu eraLjri99^ . 


Long base interferometry of fa inter y-Cas-analog s with 
V magnitudes of ss 6.6 or more (ISmith et al.ll201^ using 
the most advanced instruments currently available (e.g. 
CHARA and the ESO-VLTI) may allow us to measure 
the inclination of their circumstellar disk with enough 
accuracy so as to constrain their equatorial rotation rate 
and test the validity of this scenario. Conversely, a sys¬ 
tematic X-ray survey of the fastest rotating early Be stars 
may help to precise the critical rotation rate and effective 
temperature above which the 7 -Cas phenomenon starts 
and constrain the mechanism giving rise to the putative 
magnetic field in 7 -Cas-analogs. 


has made use of MAXI data provided by RIKEN, JAXA 
and the MAXI team. 

Facilities: RXTE (ASM), MAXI, APT at Fairborn 
Observatory. 
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XMM-Newton detection parameters and spec- 
tra of HD 120324 were extracted from the XCat-DB 
(http://xcatdb.unistra.fr/3xmm/). The observed XMM-Newton 
X-ray Hux yields a revised Log(Lx) (0.1—lOkeV) of 29.9. The 
EPIC pn X-ray spectrum is soft and requires at least three thermal 
components with a maximum temperature of Pi 2 keV. 
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APPENDIX 

A. PERIODICITIES INTRINSIC TO THE ASM 
AND MAXI DETECTORS 

In this appendix, we discuss the reality of some of the 
signals delivered by ASM/RXTE and MAXI at periods 
close to those intrinsically displayed by 7 -Cas, namely, « 
200 d (orbital period) and ~ 70d (optical quasi periodic 
variations). All timing analysis was performed using the 
least square power spectrum (LSPS) (iLomblflO^iScar^ 
Il982fl which is well suited to the unevenly distributed 
optical data while providing good results on the more 
regularly sequenced X-ray time series. 


A.I. A ^ 200 d ASM modulation artifact 

The LSPS of the entire ASM 1996-2010 broad band 
time series shows a group of peaks around 200 d tan- 
talizingly close to the 203.6 d orbital period of 7 -Cas. 
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ASM power spectra in the time intervals 1996-1999 and 
2006-2010 (dates inclusive) exhibit clear modulations at 
f»193d, while no such peak is present in between. A 
LSPS applied to the merging of the two time intervals 
listed above yields a main peak at 192 d with a weaker 
201 d alias (see Fig. [ 8 ]). With a probability of 4 x 10“® 
to be detected at random over the 20 d to 1000 d range 
of periods searched, the 192 d signal is highly signif¬ 
icant with full amplitudes in the range of 0.10 to 0.14 
cts/s. Unfortunately, this periodicity is also close to half 
a year and should therefore be considered with caution. 
Although the window spectrum does not show evidence 
of a signal at these periods, the background level does 
show a clear annual variation. Taken separately, both 
negative and positive rate time series display consider¬ 
able power at a period of exactly six months (182.6d). 
It seems thus clear that the actual background used to 
correct the light-curves varies with a half year period 
and that this modulation could impact our ability to de¬ 
tect a genuine astrophysical signal close to the 203.6 d 
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Figure 8. The Lomb-Scargle power spectrum all ASM/RXTE 
broad band data grouped in 10 d bins of time. The two higher 
peaks are at 192 d and 201 d. The dashed line shows the position 
of the orbital period of 7 -Cas. 



Figure 9. The Lomb-Scargle power spectrum of 7 -Cas MAXI 
broad band data grouped in 5 d bins of time. The highest peak 
is at 70.5 d, a period close to the precession period 70 d) of 
the orbit of the International Space Station on which the MAXI 
instrument is installed. 


orbital period of 7 -Cas. However, the reason that the 
detected period is slightly longer than 6 months is un¬ 
clear and may be due to a beating phenomenon between 
the half year per iod and a longer ^ 10 yr modulation. 
iWen et ^ (1200611 have investigated in great details the 
overall periodicities found in ASM light-curves and iden¬ 
tified a number of likely spurious frequencies. In partic¬ 
ular, periods in the range of 163 d to 196 d are detected 
in 17 sources and were therefore considered as likely ar- 
tifact s. Using power spe ctra of AGN or supernova rem¬ 
nant, [LeymeleOI] (1201 If) show that spurious signals may 
be found in the period range of 174.5 d to 191.2 d. There¬ 
fore, the occurrence of similar periodicities in the X-ray 
light curves of other ASM sources and the strong back¬ 
ground modulations with a half-year period, all point at 
a non astrophysical origin of 192 d modulation and of 
its 201 d alias. Consequently in this paper, we use ASM 
light-curves corrected for the '^200d periods. However, 
we systematically checked that using original data with¬ 
out corrections for the spurious 200 d period yield results 
consistent with those obtained based on corrected time 
series. 

A.2. A 70 d MAXI modulation artifact 


The power spectrum of the MAXI 2—20keV light 
curve of 7 -Cas (Fig. exhibits a relatively strong peak 
at P = 70.5 d, close to the « 70 d periodicity ofte n present 
in the optical photometry (jSmith et al.ll2006ll . Unfor¬ 
tunately, the best MAXI period is also identical to the 
70 d precession perio d of the International Space Station 
(|Tomida et al.l[ 2 OIIII and this obviously casts doubt on 
the reality of this detection. We note in passing that no 
evidence of a 203.6 d 7 -Cas orbital period is present in the 
power spectrum of the MAXI data, thus confirming the 
spurious nature of the signal detected in the ASM data. 
In order to investigate to which extent the precession 
period of the ISS was imprinted in MAXI light curves, 
we extracted from the MAXI database light curves of 
the 43 X-ray sources located within 40° of 7 -Cas, adding 



Figure 10. Average Least Square Power Spectrum of the 11 
MAXI sources with continuous light curves and with mean count 
rates similar to that of 7 -Cas (0.01-0.05 cnt/s). This result can be 
compared to Fig. |9] . 


another 8 sources located within 60 degrees and finally 
3 bright low-mass X-ray binaries at larger angular dis¬ 
tances from 7 -Cas (2S 0921-63, GX5-1 and LMCX-2). 
In all cases, we downloaded data averaged over one ISS 
orbital period and processed them applying the same MI¬ 
DAS procedures as used for 7 -Cas. In particular, we re¬ 
binned all individual orbit 2—20keV measurements in 
5 d bins and computed the LSPS of all light curves. We 
systematically discarded data obtained after MJD 56400 
as those acquired at a later date seem affected by stronger 
errors. About half of the extracted light curves display 
periodic gaps that appear every 70 d when the source 
moves away from the MAXI field of view. We found that 
most, if not all light curves are significantly modulated at 
the ISS precession period. Selecting the highest power in 
the period range of 66 d to 78 d, we computed the (half) 
amplitude of the best fitting sine and the mean count 
rate for all sources. Figs [TT] a & b show that the rela- 
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live (half) amplitude of the modulation does not seem 
to vary with the angular distance to 7 -Cas. However, as 
shown on Fig. [TT] c, the relative amplitude steeply de¬ 
creases with increasing source mean count rate. Such a 
behavior suggests that a background varying with the 
ISS precession period is at the origin of the spurious 
modulation. Light curves with gaps occurring every 70 d 
display systematically higher relative amplitudes due to 
their peculiar periodic window function. 7 -Cas remains 
continuously visible to MAXI at least until MJD 56400. 
The amplitude of the 70.5 d modulation of 7 -Cas is well 
within the range displayed by sources with similar X-ray 
brightness (see Fig. [TT]c). In addition, the « 70d peak 
in the average LSPS of the II MAXI sources with fluxes, 
visibility and location on the sky similar to those of 7 - 
Cas (Fig. [T0)l is virtually the same as the 70.5 d peak 
of 7 -Cas itself (Fig. |9l). We therefore conclude that the 
70.5 d modulation seen in the MAXI X-ray light curve of 
7 -Cas is mostly an artifact. We note, however, that close 
inspection of the optical and X-ray light curves indicates 
that X-ray power is maximum at times when strongest 
optical modulations occur so that a real correlation may 
exist blended with the underlying modulation produced 
by the precession of the ISS. This is supported by the 
strong close to zero time delay correlation between the 
optical and MAXI X-ray light curves observed early 2012. 


Figure 11. Dependency of the relative amplitude of the 70 d pe¬ 
riod modulation with Right Ascension (panel a), Declination (panel 
b) and mean count rate (panel c) for our sample of 55 MAXI 
sources. Count rates are in units of ph/s/cm^. Filled symbols: 
light curves without gaps. Open symbols: light curves showing 
period gaps. The position of 7 -Cas is marked by a large filled red 
lozenge. 





